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Abstract 

The production of mesons containing strange quarks (Kg, 0) and both singly and doubly strange 
baryons (A, A, and + E^) are measured at central rapidity in pp collisions at ^/s = 0.9 TeV 
with the ALICE experiment at the LHC. The results are obtained from the analysis of about 250 k 
minimum bias events recorded in 2009. Measurements of yields (dN/dy) and transverse momen- 
tum spectra at central rapidities for inelastic pp collisions are presented. For mesons, we report 
yields {{AN/dy)) of 0.184 ± 0.002 {stat.) ±0.006 [syst.) for and 0.021 ±0.004 [stat.) ± 
0.003 [syst.) for 0. For baryons, we find {AN/dy) =0.048 ±0.001 [stat.) ±0.004 {syst.) for 
A, 0.047 ±0.002 (ifaf.) ±Q.QQ5 {syst.) for A and 0.0101 ± 0.0020 (^faf.) ± 0.0009 (ijif.) for 
+ E^. The results are also compared with predictions for identified particle spectra from QCD- 
inspired models and provide a baseline for comparisons with both future pp measurements at higher 
energies and heavy-ion collisions. 




*See Appendix lAlfor the list of collaboration members 
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1 Introduction 

The production of hadrons at high transverse momenta in high energy proton-proton collisions is reason- 
ably well described by perturbative Quantum Chromodynamics (pQCD) in terms of hard parton-parton 
scattering (large momentum transfers) followed by fragmentation [1, 2|. However, the low-momentum 
region, where most particles are produced and which therefore contributes most to the underlying event, 
is dominated by soft interactions. In the soft regime, it has been found that particle production can be 
described effectively by models based on emission from an equilibrated system at a specific temperature 
and baryo-chemical potential, with additional accounting of conserved quantities [3, 4, 5]. It can also 
be treated in the framework of QCD inspired phenomenological models, that include multi-parton pro- 
cesses, extrapolated to very low-momentum transfers |6|. The contribution and evolution of multi-parton 
processes as a function of y/s is difficult to establish. Measurements of identified particles at the beam in- 
jection energy of the LHC and in the low transverse momentum (/jt) region, along with their comparison 
with QCD-inspired models, constitute a baseline for comparisons with higher centre-of-mass energies. 
The low /jt cutoff achievable through the low material budget, low central barrel magnetic field (0.5 T) 
and excellent particle identification (PID) of the ALICE detectors, allows an accurate measurement of 
the low momentum region at mid-rapidity. 

The differential transverse momentum yields (/jt spectra) and integrated yields at central rapidity of K^, 
0, A, A and + have been measured by the ALICE experiment during the commissioning phase 
of the LHC (December 2009) Q with the very first proton-proton collisions HI and are reported in 
this article. A sample of 250 k minimum bias pp collisions at = 0.9 TeV has been selected with 
triggers combining several fast detectors |9|. Measurements are performed using the tracking devices 
and the main PID detectors of ALICE in the central rapidity region {\y\ < 0.8). A comparison of the 
transverse momentum shapes (mass dependence and mean transverse momentum) with PYTHIA lU and 
PHOJET la is provided. 

This article is organized as follows. Section 2 presents the experimental conditions, the minimum bias 
event selection as well as a brief description of the main detectors and the associated event reconstruction 
tools used for the analysis. Section 3 is dedicated to the data analysis, including track and topological 
selections, signal extraction methods and the corresponding efficiency corrections. The determination 
of the systematic uncertainties are also described in this section. In section 4, the /jt spectra and the 
integrated yields of the studied particle species are given and compared with previous measurements and 
model predictions. Conclusions are given in section 5. 

2 Experimental set-up and data collection 

A detailed description of the ALICE experimental setup and its detector subsystems can be found in flOl. 
2.1 Main detectors and reconstruction techniques used for the analyses 

The central barrel of ALICE covers polar angles from 45° to 135° over the full azimuth. It is embedded in 
the large L3 solenoidal magnet providing a nominal magnetic field B of 0.5 T. Within the barrel, the two 
tracking detectors used in these present analyses consist of an Inner Tracking System (ITS), composed 
of 6 cylindrical layers of high-resolution Silicon detectors and a cylindrical Time Projection Chamber 
(TPC). PID is performed using secondary (displaced) vertex reconstruction, invariant mass analysis and 
single track PID methods, which include the measurement of specific ionization in the ITS and the TPC, 
and the information from the Time-Of-Flight detector (TOE). 

2.1.1 The Inner Tracking System 

The Silicon Pixel Detector (SPD) corresponds to the two innermost ITS layers. These two layers have 
a very high granularity with a total of about 9.8 million pixels, each with a size of 50 X 425 ixm^. They 
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are located at radii of 3.9 and 7.6 cm and the pseudorapidity coverages are |t]| < 2.0 and |t]| < 1.4 
respectively. The detector provides a position resolution of 12 /xm in the r<p direction and about 100 /^m 
in the direction along the beam axis. It can also deliver a signal for the first level of trigger (LO) in less 
than 850 ns. The two layers of the Silicon Drift Detector (SDD), located at radii of 15.0 and 23.9 cm, 
are composed of 260 sensors, including 133 000 collection anodes with a pitch of 294 /xm. They provide 
a charge deposit measurement and a position measurement with a resolution of about 35 /^m in the 
direction and about 25 ;Umin the beam direction lITTl . The Silicon Strip Detector (SSD) consists of 1698 
double-sided sensors (with a strip pitch of 95 /im and a stereo angle of 35 mrad) arranged in 2 layers 
located at radii of 38 and 43 cm. It provides a measurement of the charge deposited in each of the 2.6 
million strips, as well as a position measurement with a resolution of 20 /im in the T(j) direction and about 
800 umin the beam direction. 

The ITS sensor modules were aligned using survey information and tracks from cosmic-ray muons and 
pp collisions. The corresponding methods are described in ITTll . 

The percentage of operational channels in the ITS during the 2009 run is 82% for the SPD, 91% for the 
SDD and 90% for the SSD. 



2.1.2 The Time Projection Chamber 

The ALICE TPC is a cylindrical drift detector with a pseudorapidity coverage of |t]| < 0.9 (TT\. It has 
a field cage filled with 90 m^ of Ne/C02/N2 (85.7/9.5/4.8%). The inner and outer radii of the active 
volume aie of 85 cm and 247 cm respectively and the length along the beam direction is 500 cm. Inside 
the field cage, ionization electrons produced when charged particles traverse the active volume on either 
side of the central electrode (a high voltage membrane at —100 kV) migrate to the end plates in less than 
94 us. A total of 72 multi-wire proportional chambers, with cathode pad readout, instrument the two 
end plates of the TPC which are segmented in 18 sectors and amount to a total of 557, 568 readout pads. 
The ALICE TPC ReadOut (ALTRO) chip, employing a 10 bit ADC at 10 MHz samphng rate and digital 
filtering circuits, allows for precise position and linear energy loss measurements with a gas gain of the 
order of 10^. 

The position resolution in the direction varies from 1100 fim to 800 fim when going from the inner 
to the outer radius whereas the resolution along the beam axis ranges between 1250 /im and 1 100 /xm. 



2.1.3 The Time-Of-Flight detector 

The ALICE Time-Of-Flight detector fVS\ is a cylindrical assembly of Multi-gap Resistive Plate Cham- 
bers (MRPC) with an inner radius of 370 cm and an outer radius of 399 cm, a pseudorapidity range 
|rj| < 0.9 and full azimuth angle, except for the region 260 < <p < 320 at rj near zero where no TOE 
modules were installed to reduce the material in front of the Photon Spectrometer. The basic unit of the 
TOE system is a 10-gap double-stack MRPC strip 122 cm long and 13 cm wide, with an active area of 
120 X 7.4 cm^ subdivided into two rows of 48 pads of 3.5 x 2.5 cm^. Eive modules of three different 
types are needed to cover the full cylinder along the z direction. All modules have the same structure 
and width (128 cm) but differ in length. The overall TOE barrel length is 741 cm (active region). It 
has 152,928 readout channels and an average thickness of 25 — 30% of a radiation length, depending on 
the detector zone. Eor pp collisions, such a segmentation leads to an occupancy smaller than 0.02 %. 
Its front-end electronics is designed to comply with the basic characteristics of a MRPC detector, i.e. 
very fast differential signals from the anode and cathode readout. Test beam results demonstrated a time 
resolution below 50 ps, dominated by the jitter in the electronic readout. 
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2.1.4 The VZERO Counters 

The VZERO counters are two scintillator hodoscopes located along the beam direction at —0.9 m and 
3.3 m from the geometrical centre of the experiment. They correspond to a coverage of— 3.7<T7<— 1.7 
and 2.8 < tj < 5.1 respectively and have a time resolution close to 0.5 ns. They are used as trigger 
detectors and help to remove beam-gas interaction background. 

2.1.5 Track reconstruction and particle identification 

The global tracking system in the ALICE central barrel (combining the ITS and the TPC) covers the 
pseudorapidity window |?7| < 0.9. 

The reconstruction in the tracking detectors begins with charge cluster finding. The two coordinates of 
the crossing points (space points) between tracks and detector sensitive elements (pad rows in the TPC, 
and silicon sensors in the ITS) are calculated as the centres of gravity of the clusters. The errors on the 
space point positions are parametrized as functions of the cluster size and of the deposited charge. In the 
TPC, these errors are further corrected during the tracking, using the crossing angles between tracks and 
the pad rows. 

The space points reconstructed at the two innermost ITS layers (pixel detector, SPD) are then used for the 
reconstruction of the primary vertex. One space point from the first SPD layer and one from the second 
layer are combined into pairs called "tracklets". The primary vertex is consequently reconstructed in 3D 
as the location that minimizes the sum of the squared distances to all the tracklet extrapolations. If this 
fails, the algorithm instead reconstructs the z coordinate of the vertex by correlating the z coordinates 
of the SPD space points, while for x and y the average position of the beam in the transverse plane 
(measured basis by a dedicated calibration procedure on a run-by-run basis) is assumed. 

Track reconstruction in ALICE is based on the Kalman filter approach and is discussed in detail in fT5l. 
The initial approximations for the track parameters (the "seeds") are constructed using pairs of space 
points taken at two outer TPC pad rows separated by a few pad rows and the primary vertex. The primary 
vertex position errors for this procedure are considered to be as big as 3 cm. The seeds for the secondary 
tracks are created without using the primary vertex, since such a constraint would unnecessarily reduce 
the VO finding efficiency. The space points are searched along the line connecting the pairs of points 
taken at those two outer TPC pad rows. 

Once the track seeds are created, they are sorted according to the estimate of their transverse momentum 
(;7t)- Then they are extended from one pad row to another in the TPC and from one layer to another 
in the ITS towards the primary vertex. Every time a space point is found within a prolongation path 
defined by the current estimate of the covariance matrix, the track parameters and the covariance matrix 
are updated using the Kalman filter. For each tracking step, the estimates of the track parameters and 
the covariance matrix are also corrected for the mean energy loss and Coulomb multiple scattering in 
the traversed material. The decision on the particle mass to be used for these corrections is based on the 
dE/dx information given by the TPC when available. If the information is missing or not conclusive, a 
pion mass is assumed. Only five particle hypotheses are considered: electrons, muons, pions, kaons and 
protons. 

All the tracks are then propagated outwards, through the ITS and the TPC. When possible, they are 
matched with the hits reconstructed in the TOE detector. During this tracking phase, the track length 
and five time-of-flight hypotheses per track (corresponding to the electron, muon, pion, kaon and proton 
masses) are calculated. This information is later used for the TOE PID procedure. The track parame- 
ters are then re-estimated at the distance of closest approach (DCA) to the primary vertex applying the 
Kalman filter to the space points already attached. Einally, the primary vertex is fitted once again, now 
using reconstructed tracks and the information about the average position and spread of the beam-beam 
interaction region estimated for this run. 
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In pp collisions, the track reconstruction efficiency in the acceptance of TPC saturates at about 90% 
because of the effect of the dead zones between its sectors. It goes down to about 75% around pj = 
1 GeV/c and drops to 45% at 0.15 GeV/c. It is limited by particle decays (for kaons), track bending at 
low Pi and absorption in the detector material. The amount of material traversed by particles near rj =0 
is about 11% of of a radiation length including the beam pipe, the ITS and the TPC (with service and 
support). 

The overall pj resolution is at least as good as the TPC-standalone resolution, which is typically 1 % for 
momenta of 1 GeV/c and 7% for momenta of 10 GeV/c, and follows the parameterization (cj(pt) / Pt)^ = 
(0.01)^ + (0.007 pjy where pj is expressed in GeV/c (see llT4ll for the details). 

The resolution of the track transverse impact parameter (the minimal distance between a track and the pri- 
mary vertex in the transverse plane) depends on the precision of track and primary vertex reconstruction. 
These in turn depend on the momentum, and, in the case of the vertex, on the number of contributing 
tracks. As it was estimated from the data, the transverse impact parameter resolution for a typical pp 
event could be parameterized as a{pj) = 50 + 60/ (/^t)"'^ (c is in /xm, and pj is in GeV/c), which was 
defined by the level of the ITS alignment achieved in 2009. 

The dE/dx resolution of the TPC is estimated to be about 5% for tracks with 159 clusters lfT2l . which is 
better than the design value |[T5l . When averaged over all reconstructed tracks, this resolution is about 
6.5%. 

During the run, the preliminary calibration of the TOP detector corresponds to a resolution of 180 ps, 
which includes 140 ps due to the jitter in the absolute time of the collisions. This contribution is reduced 
to about 85 ps for those events with at least 3 tracks reaching the TOP, in which case an independent 
time zero determination is possible. The matching efficiency with TPC tracks (which includes geometry, 
decays and interaction with material) is on average 60% for protons and pions and reaches 65% above 
Pi = I GeV/c. For kaons it remains sligthly lower [16|. Above pi = 0.5 GeV/c, the TOP PID has an 
efficiency larger than 60% with a very small contamination. 

2.2 LHC running conditions and triggers 

For the first collisions provided by the Large Hadron Collider, four low intensity proton bunches (10^ 
protons per bunch, giving the luminosity of the order of 10^^ cm^^s^^) per beam were circulated, and two 
pairs of them crossed at the ALICE interaction point. Under such conditions, the rate for multiple events 
in a given bunch-crossing ("pile-up") was negligible. The energy in the centre of mass corresponded to 
twice the beam injection energy, that is = 0.9 TeV. The data acquisition of ALICE was triggered by 
requiring two coincidence conditions: i) the LHC bunch-crossing signal together with the two beam pick- 
up monitors (BPTX); ii) ALICE minimum bias (MB) trigger requiring a combination of signals from the 
SPD and from the VZERO counters. For these analyses, the MBqr was used, which is fulfilled when at 
least one of the VZEROs or the SPD trigger is fired [91. The corresponding data rate was ~ 10 Hz. 

3 Data analysis 

3.1 Event and track selection 

The primary vertex is reconstructed using either SPD tracklets IH (5% of the events) or global tracks 
(95% of the events). Events are selected by requiring that the distance between the position of primary 
vertex and the geometrical centre of the apparatus along the beam axis be less than 10 cm (z = —0.40 cm 
and rmSj = 4.24 cm, where z is the average position of the primary vertex along the beam axis). Events 
with less centred primary vertices (|z| > 10 cm) are discarded in order to minimize acceptance and 
efficiency biases for tracks at the edge of the TPC detection volume. The average position and dispersion 
for both horizontal and vertical directions are found to be I = —0.35 mm (y = +1.63 mm) and rmsj^ = 
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Fig. 1: Primary vertex distributions for the analysed events. The left panel shows the distributions along the beam 
axis. Selected events (full symbols) are required to have a reconstructed primary vertex with |z| < 10 cm. The right 
panel corresponds to the directions perpendicular to the beam axis: horizontally (i.e. x-direction, squares and full 
line) and vertically (i.e. y-direction, triangles and dashed line). 



0.23 mm (rms^ = 0.27 rmn). No conditions were applied on the x and y position of the vertex. The 
total number of events used for obtaining the particle spectra and yields is about 250 k events. Figure [T] 
shows the primary vertex distribution along the beam axis (left panel) and for the x and y directions (right 
panel). The dashed lines indicate the limits of the selected vertex region. 

The normalization to the number of inelastic events (INEL) is obtained in the same way as other ALICE 
analyses HI [161. H leads to a correction for the normalization of ~ 5% with an uncertainty of 2%. This 
uncertainty is added to the ones described in section 1331 and mainly related to the modeling of the fraction 
of diffractive events with several Monte Carlo event generators. 

Several quality criteria are defined for track selection. Each track is required to have been reconstructed 
in the TPC in the initial outward-in step of tracking and then successfully refitted in the final back- 
propagation to the primary vertex as described in section 12.1.51 It is also required that each track has at 
least 80 TPC clusters out of a maximum of 159. At the reconstruction level, split tracks are rejected as 
well as those which may correspond to daughters of kaons decaying in the TPC. 

As the particle is a strongly decaying resonance, its daughters are indistinguishable from primary 
particles at the reconstruction level and therefore primary track selections are used. As a first step, each 
track is propagated to the reconstructed primary vertex. If this operation is successful, the track is kept 
if it has a DCA smaller than 5 mm (3 cm) in the transverse (longitudinal) direction with the additional 
constraints of having at least one SPD cluster and a of less than 4 per cluster assignment (for each 
cluster, the has two degrees of freedom). 

Depending on its hfetime, a particle may cross several layers of the ITS before weakly decaying. The 
probability that the daughter tracks of Kg, A, A and + have a hit in this detector decreases accord- 
ingly. Therefore, no specific condition on the number of ITS hits is required for the daughter tracks of 
the reconstructed secondary vertices. However, other quality criteria are applied for selecting the daugh- 
ter tracks of weakly decaying particles which are not considered as primaries. The selections described 
are summarized in Table [T] The measurement of differential yields in rapidity and pi bins cannot be 
performed simultaneously for the particles considered due to the small available statistics. Therefore 
the rapidity ranges are chosen such that i) the efficiency does not vary strongly for each species and ii) 
the rapidity distribution is sufficiently flat for it to be possible to rely on the Monte Carlo to obtain the 
corrections. 
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Table 1: Track selection criteria. 



Common selections 


Detectors required for track rec/fit 




ITS,TPC 


Number of TPC clusters ^ 


> 


80 


N((7) dfi/dx (TPCPID) 




2, to 5 


Primary track selections 


per cluster 


< 


4 


DCA to primary vertex (r,z) 


< 


(0.5,3.0) cm 


Number of SPD clusters 


> 


1 


Secondary track selections 


transverse momentum 


> 


160 MeV/c 


^ maximum number for the TPC is 159; 






maximum number for the SPD is 2; 






in the cases of Kg, A and A. 







Table 2: Secondary vertex selection criteria. 



Common selections 



Minimum transverse decay radius 


> 


0.2 cm 


Maximum transverse decay radius 


< 


100 cm 


yO vertex selections (Kg, A and A) 


DCA of VO daugther track 






to primary vertex 


> 


0.05 cm 


DCA between VO daughter tracks 


< 


0.50 cm 


Cosine of VO pointing angle (A and A) 


> 


0.99 


Cascade vertex selections 


DCA of cascade daughter track 






to primary vertex " 


> 


0.01 cm 


DCA between VO daughter tracks 


< 


3.0 cm 


Cosine of VO pointing angle 


> 


0.97 


DCA of VO to primary vertex 


> 


0.001 cm 


VO invariant mass 


> 


lllOMeV/c^ 


VO invariant mass 


< 


1122MeV/c2 


DCA between VO and bachelor track 


< 


3.0 cm 


Cosine of cascade pointing angle 


> 


0.85 



for bachelor and each VO daughter. 
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Table 3: Main characteristics of the reconstructed particles: valence quark content, mass, cT and charged decay 
branching ratio (B.R.) IJ8J . 





Particles 


mass (MeV/c^) 


CT 


charged decay 


B.R. (%) 


Mesons 




497.61 


2.68 cm 




69.2 


{ss) 


1019.46 


45 fm 




49.2 


Baryons 


A {uds) and A (uds) 
(dss) and {dss) 


1115.68 
1321.71 


7.89 cm 
4.91 cm 


A ^ p + K^ and A — > p + 7r+ 
E" ^ A + TT^ and — > A + 7r+ 


63.9 
99.9 



3.2 Particle reconstruction and identification methods 
3.2.1 Topological reconstruction o/K^, A, A and S + 

TheK'^, A, AandE^ + S are identified by applying selections on the characteristics of their daughter 
tracks (see Table O and using their weak decay topologies in the channels listed in Tabled 

The measurement of K^, A and A is based on the reconstruction of the secondary vertex (VO) associated 
to their weak decay. The VO finding procedure starts with the selection of secondary tracks, i.e. tracks 
having a sufficiently large impact parameter with respect to the primary vertex. All possible combinations 
between two secondary tracks of opposite charge are then examined. They are accepted as VO candidates 
only if the DCA between them is smaller than 0.5 cm. The minimization of the distance between the 
tracks is performed numerically using helix parametrizations in 3D. The VO vertex position is a point 
on the line connecting the points of closest approach between the two tracks. Its distance from each 
daughter track is taken to be proportional to the precision of the track parameter estimations. Once their 
position is determined, only the VO candidates located inside a given fiducial volume are kept. The 
inner boundary of this fiducial volume is at a radius of 0.2 cm from the primary vertex, while the outer 
limit is set at 100 cm. Finally, for A and A reconstruction, the VO finding procedure checks whether the 
particle momentum {p) associated with the VO candidate (calculated as the sum of the track momenta 
extrapolated to the position of the DCA) points back to the primary vertex. This is achieved by applying 
a cut on the cosine of the angle (pointing angle 6p) between p and a vector connecting the primary vertex 
and the VO position (cos Qp > 0.99). The invariant mass of each candidate can then be calculated either 
under the or the A hypothesis. 

The TPC FID helps substantially to remove the combinatorial background for the A and A (mainly for 
the baryon daughter identification, while it is not needed for the decaying into pions). TFC FID 
is described in paragraph 13.2.31 The selections here concern the proton daughter only and have been 
chosen to be looser for the daughter track with momentum below 0.7 GeV/c (±5a) and tighter for 
higher momentum (ib3a). 

The + E^ particles are identified via their "cascade" decay topology. The cascade finding procedure 
starts from the VO finding procedure for the A daughter but with less stringent selection criteria (see 
Table |2] and Cascade vertex selections). This is done to increase the efficiency and to allow for the fact 
that the daughter A's do not have to point back to the primary vertex. 

The VO candidates found within the A mass window (1116 ±6 MeV/c^) are combined with all pos- 
sible secondary tracks (bachelor candidates) with the exception of both VO daughter tracks. A cut on 
the impact parameter of the bachelor track is applied to reject the primary particles which increase the 
combinatorial background. 

A VO-bachelor association is performed if the distance of closest approach between the bachelor track 
and the VO trajectory (DCA between VO and bachelor track) is small (less than 3 cm). Finally, this 
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Fig. 2: Invariant mass distributions of Kg, A and A, (j) and the sum Z +Z . The vertical arrows indicate the 
nominal mass values from PDG. 

cascade candidate is selected if its reconstructed momentum points back to the primary vertex (cosine of 
cascade pointing angle). The cascade finding is limited to the fiducial region used for VO reconstruction 
(see Table 121). 

In addition to topological selections, the reconstruction of cascades also makes use of the single-track 
PID information delivered by the TPC. This is considered for each of the three daughters (both pions 
and the proton). For each track, a loose selection is required (ib4a over the whole momentum range) to 
reject the combinatorial background in part. The resulting invariant mass distributions are presented in 
Fig. El 



3.2.2 Additional quality checks for Kg, A, A 

A significant fraction of the reconstructed VO come from 7 conversion in the detector material. This can 
be clearly seen in the Armenteros-Podolanski distribution |17 | shown in Fig. [3l where and are 
the longitudinal components of the total momentum for the positive and negative daughters respectively, 
relative to the direction of the VO momentum vector. The Kg, A and A signal regions are symmetric and 
clearly distinguishable. 

The lifetime (ct) distributions for K^, A and A are also checked. All VO candidates within a ib3a 
effective mass region around the nominal value are used in the distribution without further residual 
background subtraction. The corresponding distributions of ct = are obtained, where L is defined as 
the distance between primary and VO vertices, and m and p are the particle mass and momentum. Because 
of the acceptance, the single track efficiency and the topological selections applied at reconstruction 
level, the reconstruction efficiency as a function of the decay length is not constant. The corresponding 
corrections are extracted from the reconstruction of full Monte Carlo simulations (see section 1X41 ). The 
corrected cT distributions are fitted using exponential functions. The results are shown with the statistical 
uncertainties in Fig.lH The extracted decay lengths of 7.9 ±0.1 cm, 7.7 ±0.1 cm and 2.72 ±0.03 cm for 
A, A and Kg, respectively, are compatible with the PDG values given in Tabled 
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Fig. 3: Armenteros-Podolanski distribution for VO candidates showing a clear separation between Kg, A and A. 
The 7 converting to e+e^ with the detector material are located in the low qr region, where qr is the momentum 
component perpendicular to the parent momentum vector. 




Lffi(cm) 

Fig. 4: Kg, A and A hfetime distributions obtained for the candidates selected by the invariant mass within a ±3(7 
region around the nominal mass and corrected for detection efficiency. The distributions are scaled for visibihty 
and fitted to an exponential distribution (straight fines). Only statistical uncertainties are shown. 
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momentum p (GeV/c) momentum p (GeV/c) 

Fig. 5: Left panel shows the distribution of the measured energy-deposit of charged particles versus their momen- 
tum in the TPC. The dashed lines delimit successively a ±5(7 then a ±3(7 selection of the kaon tracks using the 
ALEPH parameterization lfT9l of the Bethe-Bloch curve (solid line). Right panel shows the relative difference 
between TOP measured times and that corresponding to a kaon mass hypothesis. The dashed lines delimit a coarse 
fiducial region compatible with this kaon hypothesis. 

3.2.3 ^ reconstruction 

The resonance is reconstructed through its principal decay channel <p — K+K^ (see Table [3]). With 
a cT of 45 fm, its decay vertex is indistinguishable from the primary collision vertex. Therefore the 
selection criteria adopted for the candidate daughter tracks are the ones used for primaries, as specified 
in Tabled 

A crucial issue for the <p reconstruction, as for any strongly decaying resonance, is the combinatorial 
background determination. In the present analysis PID is used to select kaons, rejecting most of the 
background while leading to a very small loss in efficiency. For this purpose, tracks are selected if the 
PID information from the TPC is compatible with a kaon signal and using the TOP signal when available. 

For each track, the expected energy loss is calculated using a parametrised response based on the Bethe- 
Bloch formula fT9l computed with a kaon mass hypothesis. It is compared with the TPC specific ion- 
ization dE/dx measured via truncated mean (the reconstructed momentum being evaluated at the inner 
radius of the TPC). With the current TPC calibration for this data set, the assumed dE/dx resolution is 
6 %. For momenta smaller than 350 MeV/c, the species are well separated so the window is set to ±5a 
with little or no contamination; above 350 MeV/c, it is set instead to ±3a as shown in the left panel of 
Fig.S 

The accepted band for TOF kaon identification is defined with two hyperbolas as shown in the right panel 
of Fig. [51 

3.3 Background evaluation and signal extraction 

For minimum bias pp collisions, the signals for all particles are clearly distinguishable from the combina- 
torial background as shown in Fig. [2] Two different methods are used to extract the invariant mass signal 
from the background. For the single strange particles (K*^, A and A), the signal is first approximated by 
a Gaussian on a second order polynomial background. This gives an estimate of the signal mean and 
width although the invariant mass signal is not strictly Gaussian. Then the background is sampled on 
each side of the signal by using both sampled regions that are more than 6a away from the Gaussian 
mean. The assumption that no reconstructed signal is included in these regions is checked using Monte 
Carlo data. The width of the background regions can vary depending on the pj interval considered in the 
invariant mass distributions. The sum of signal and background (S+B) is sampled in the region defined 
by the mean ±4a. 
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Fig. 6: Plot illustrating the "bin-counting" method used to extract the raw yields. It corresponds to the invariant 
mass distribution of Kg for the pj bin [0.4 — 0.5] GeV/c. The hashed regions show where the background is 
sampled; they are chosen to be 6(7 away from the signal approximated with a Gaussian distribution. The averaged 
or fitted background is subtracted from the signal region of ±4a. 

The sampling method is illustrated in Fig.[6]for the K^. Two methods are used to evaluate the background 
and give consistent results. The background areas are either i) fitted simultaneously with polynomial 
functions (from first to third order) or ii) averaged by simply counting the number of entries ("bin- 
counting"). The background B under the signal S is estimated using the normalized area sampled on 
both sides of the signal region (Gaussian mean ib4a). The signal yield S = (B + 5) — B is thus evaluated 
without any assumption as to its shape. Systematic effects such as signal asymmetry are taken into 
account by varying the size of the signal and background intervals up to la. The difference between the 
two methods (fit and bin-counting) contributes to the evaluation of the systematic uncertainties associated 
to the signal extraction. 

In the case of the + S^, statistical uncertainties are significant so that, in parallel to the bin-counting 
method, the background level is simply estimated by a straight line fit. 

The invariant mass distribution has a larger combinatorial background and a function reproducing both 
the background and the signal is preferred. It is found that the background can be well reproduced by 
a function f(M) = a\/M — b, while the peak has the shape of a Gaussian. The peak range is defined as 
ib4a around the PDG mass of the 0, where a = r/2.35 and Y is the nominal value of the resonance full 
width at half maximum (4.5 MeV/c^) [ 18 |. For each analyzed pi bin, several fit ranges are investigated. 
It is found that the fitted width matches that extracted from a full Monte Carlo simulation (as defined in 
section l34l ) within 5%, except for the last p-^ bin where it is broader (~ 10%). While fluctuations of the fit 
values as a function of the fit range are taken into account for the systematic error (see section |3".5. IK the 
fit values used for all subsequent steps in the analysis are those that minimize the difference |;^^/NDF — 
1|. Figure |7] illustrates the method for the [1.0— 1.5] GeV/c pt: bin. Every unlike-sign track pair passing 
all selection criteria and falling within the invariant mass peak range is counted. The total number of 
is estimated by subtracting the integral of the background function alone, computed in the same invariant 
mass range. 



The signal counts (raw yields) for each of the pi bins are histogrammed as a function of pi for K^, A, 
A in Fig. [8] and for and + in Fig. |9l The uncertainties correspond to both the statistical errors 
related to the number of counts and the systematics from the bin-counting and fit methods used to extract 
the signal from the background. 
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Fig. 7: Background evaluation for the (j) corresponding to the pj bin [1.0-1.5] GeV/c. The inset shows the (j) signal 
after background subtraction. The vertical arrows indicate the nominal mass value from PDG. 
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Fig. 8: Reconstructed (raw) yields of Kg (open circles), A (open squares) and A (open triangles) as a function of 
Pt- The change of bin size results in successive offsets of the raw counts at pj — [1 .4, 1 .6, 2.4] GeV/c for Kg and 
pr = [1.6,2.4,3.0] GeV/c for A and A. Uncertainties correspond to the statistics and the systematics from the 
signal extraction. They are represented by the vertical error bars. The horizontal error bars give the bin width. 
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Fig. 9: Reconstructed (raw) yields of (stars) and Z +Z (diamonds) as a function of pj. With the current 
statistics, 4 pj bins ai-e used for the <j) ([0.7-1.0], [0.7-1.5], [1.5-2.0] and [2.0-3.0] GeV/c) and 3 pj bins for the 
+ ([0.6-1.4], [1.4-2.0] and [2.0-3.0] GeV/c). Uncertainties correspond to the statistics (i.e. the number of 
reconstructed particles) and the systematics from the signal extraction. They are represented by the vertical error 
bars (the horizontal ones give the bin width). 
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Fig. 10: Efficiency of (open circles), A (open squares) and A (open triangles) as a function of pi. The 
uncertainties correspond to the statistics in Monte Carlo samples used to compute the corrections. The efficiency 
is limited by the branching ratio represented by a solid arrow for (0.692) and by a dashed arrow for A and A 
(0.639). 
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Fig. 11: Efficiency of ^ (stars) and + (diamonds) as a function of pj. The uncertainties correspond to the 
statistics in the Monte Carlo sample used to compute the corrections. The efficiency is limited by the branching 
ratio represented by a solid arrow for + E (0.636) and by a dashed arrow for (0.492). 



3.4 Efficiency corrections 

The efficiency corrections are obtained by analysing Monte Carlo (MC) events in exactly the same way 
as for the real events. Little dependence is found on the several MC generators which are used. Therefore 
the corrections presented here are obtained using the event generator PYTHIA 6.4 (tune D6T) HJ |20]| 
and GEANT3 Il2ll for particle transport through the ALICE detectors. 

The MC information is propagated through the whole reconstruction and identification procedure to 
generate the differential pj efficiencies as shown in Fig. [TO] for Kg, A and A and in Fig. [TT] for and 
S~ + S^. The uncertainties correspond to the statistics of Monte Carlo samples used to compute the 
corrections. For all particles, the global efficiency is limited at low p^ because of the acceptance of 
at least two charged daughter tracks in the detection volume of the TPC (three tracks in the case of 
+ It rapidly increases with pi but cannot exceed the asymptotic limits given by the charged 
particle decay branching ratios presented in Table |3] The difference between the A and A reflects the 
absorption of the anti-proton daughter of the A. For all the variables used to select the particles and 
improve the signal over noise ratio (see Tables [J and |2l), it is verified that data and MC distributions 
match, thus possible efficiency biases can be properly managed. Examples of such distributions are 
presented in Fig. [T2j 
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Fig. 12: Comparison between data (red circles) and Monte Carlo (black open triangles) for several topological vari- 
ables used to select secondary vertices. The top panels correspond to the candidates selected in a ±20 MeV/c^ 
invariant mass window around the nominal mass. The distribution of the DCA between the positive daughter track 
and the primary vertex and the DCA distribution between the two daughters tracks are displayed in the left and the 
right top panels respectively. On the bottom panels, the same distributions are shown for the A candidates selected 
in a ±8 MeV/c^ invariant mass window around the nominal mass. 
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Table 4: Point-to-point systematic uncertainties expressed in percentage for pj spectra of different particles. For 
each particle, the reported values correspond to the effect on the lowest pj bin, the average and the highest pj bin, 
except for the feed-down contributions where values are estimated as being constant versus pi or where the effect 
is found to be negligible (less than 2 standard deviations from the default value on the corrected spectrum). 



systematic effects (%) 
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3.5 Estimation of the systematic uncertainties 

Systematic uncertainties are discussed in the following sections, where details are given on the contri- 
butions due to topological selections and signal extraction methods, as well as those due to material 
budget and feed-down. As for efficiency corrections, MC data are generated with PYTHIA 6.4 (tune 
D6T) lUlEOI and transported with GEANT3 li2ri . At low pT> the anti -proton absorption cross section 
in GEANT3 is known to be too large ||22l|23l. GEANT4 (with the absorption cross-sections of [24]) 
was then used to correct the anti-proton tracking efficiency. The information is summarized in Table |4] 
In addition to these point-to-point systematic uncertainties, there is also a 2 % systematic error on the 
global normalization coming from the evaluation of the total number of inelastic events. 

3.5.1 Systematic uncertainties due to track or topological selections and signal extraction 

Systematic uncertainties due to tracking and topological identification are determined by varying the 
track and topological (for secondary vertices) selections, as well as the definition of the regions sampled 
for signal extraction. To assess the different systematic uncertainties, only the deviations that are statisti- 
cally significant are taken into account (more than 2 standard deviations away from the central value on 
the corrected spectrum). 

The systematic variation of track and topological selections results in a variation of the amount of signal 
extracted from invariant mass distribution in both data and the Monte Carlo simulation mentioned above. 
The difference between these amounts of signal corresponds indirectly to the accuracy with which the 
MC simulation reproduces the characteristics of real events, from the simulation of the detector response 
to the background shape and composition considered for the extracted signal. It is estimated that the 
point-to-point uncertainties in the pi spectra are at most 4.6 %, 3.3 %, 4.7 %, 6 % and 13.9 % for the 
Kg, A, A, and + S^, respectively. The systematic uncertainties of the signal extraction for K^, A, A 
and + are /^T-dependent and estimated by varying the invariant mass regions where the signal and 
the background are sampled using the bin-counting method described in section 13.31 For the signal, 
the systematic uncertainties from background subtraction are estimated using three different criteria. 
First of all, the function reproducing the background is replaced by a second or third order polynomial. 
Moreover, the fit is repeated fixing the width parameter to ±10% of the value obtained in the default 
procedure (described in section [331) . and also to the value obtained when fitting the Monte Carlo sample 
and to ±10% of this. Finally, the fit range is also varied. All of these computations result in a variation 
of the raw counts with respect to those shown in Fig. |9] Although a quite large compatibiUty region is 
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requested for PID (at least 3 a) the effects of varying the dE/dx selections are taken into account for the 
corresponding efficiency calculation. For both and + S , statistical errors dominate after signal 
extraction (see section 13.31 ) and consequently, some systematic effects due to PID are extrapolated from 
single track and VO measurements. The TOF PID selection is applied only to reject the (/> background. No 
systematic effects are observed on the signal i) for the Monte Carlo data sample, when the selection is 
applied to the daughters in addition to all other cuts; ii) for real events, when comparing the statistics 
before and after applying the selection. 

3.5.2 Systematic uncertainties due to material budget and absorption cross-section 

A dedicated study involved the variation of the detector material thickness crossed by particles. The 
material budget uncertainty, based on y conversion measurements, is estimated to be 7% in terms or 
radiation length |22|. The efficiency variation due to this material budget uncertainty depends on the 
momentum of each of the decay daughters. Although such a variation is also correlated with the mo- 
mentum of the parent particle, the corresponding systematic uncertainties are reported as point-to-point 
enors in Table |4]for the lowest, the average and the highest /jt bin and eventually added in quadrature to 
the total systematic errors. 

Specific uncertainties are related to the (anti-)proton absorption and scattering cross-sections used for 
propagating these particles through the geometry of the detectors with both GEANT3 ||2T]| (and its default 
absorption cross-sections) and GEANT4 (using the absorption cross-sections of 1*241). More details about 
the modifications can be found in |[22l l23l and references therein. The corresponding corrections are 
taken into account in the efficiency versus /?t assuming that absorption cross-sections are identical for 
the (anti-) hyperon and its (anti-) proton daughter. The uncertainties associated with these corrections are 
derived from the (anti-)proton cross-section uncertainties and the values are estimated as constant and 
lower than 1% (2%) for A (A) and 2% for + 

3.5.3 Systematic uncertainties for A and A due to feed-down 

Some of the reconstructed A and A particles come from decays of S-hyperons. The proportion of recon- 
structed secondary A and A depends on the selection criteria used. For the parameters listed in Table |2] 
(VO vertex part), the impact of the E feed-down on the final spectra is evaluated to be 13% for A and 
12% for A. No pt: dependence is found within uncertainties. 

This assessment results in a global correction of the spectra, applied as an additional factor in the overall 
normalization. Provided that both primary and secondary A have similar spectral shapes, such integrated 
correction is applicable. This is tested directly using Monte Carlo data, but also with real data, changing 
the fraction of the secondary A by varying the DCA of reconstructed candidates. Within the available 
statistics and pi reach, no significant change in spectral shape is observed. 

Using Monte Carlo, the ratio rfeed-down of the reconstructed (S^) candidates to the number of recon- 
structed A (A) candidates from E decays is: 

{N'p-)mc 

''feed— down 



(^A^E-)mc 

Assuming that this ratio is the same in both Monte Carlo and data, the whole feed-down contribution to 
the spectra is estimated by dividing the number of reconstructed (S^) in data by the ratio extracted 
from Monte Carlo: 

(^A^S-)data = -^^ 

'feed— down 

Besides the E contribution, other sources may feed the A population resulting in additional systematic 
uncertainties. In Monte Carlo simulations, A particles possibly generated in the detector material induce 
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Fig. 13: Comparison of the corrected yields as a function of pj for Kg (circle) and charged kaons (K+) (open 
squares), identified via energy loss in the TPC and ITS, and via time of flight in the TOR The points are plotted at 
the centre of the bins. The full vertical lines associated to the Kg points, as well as the gray shaded areas associated 
to the K+ points, correspond to the statistical and systematic uncertainties summed in quadrature whereas the 
inner vertical lines contain only the statistical uncertainties (i.e. the number of reconstructed particles) and the 
systematics from the signal extraction. 

a 1 .7 % uncertainty. The same uncertainty in the case of A is below 1 %. The contiibution from Q. decays 
is found to be negligible. It should be noted that since A (A) from electromagnetic iP (£*') decays cannot 
be distinguished from the direct ones, the identified A (A) include these contributions. 

3.6 Pi spectra and global yield extraction 

The Kg spectrum is first shown on a linear scale in Fig. [13] and compared with charged kaon spectra fT6l. 
Within uncertainties, good agreement is found between K^ and K+ in the measured range. 

Figure Opresents the corrected p-^ spectra for all species, including both statistical errors and systematic 
uncertainties. The spectra are fitted with two different functional forms in order to extract the global 
integrated yields: 



dydpi 
d^N 
dydpi 



Ax pjX e T 

{n-\){n-2) dN 
nT[nT +m{n — 2)] dy 
ntT — m 



(1) 



XpjX 1 + 



riT 



(2) 



where mj = ym^ + pj. The pj exponential has two parameters: the normaUzation A and the inverse 
slope parameter T. The Levy function [Eq. ©J, already used at lower energies (25), is shown to be useful 
when the pj range is wide: it includes both an exponential shape for low pj (which can be characterized 
by an inverse slope parameter T) and a power law component (governed by the power parameter n) for 
the higher pj region. The results of these fits to the spectra, where statistical and systematic uncertainties 
are added in quadrature, are shown in Fig. [14] and in Table [5] In the case of K^ for which the statistics 
and the pj range are larger than for other species, the ;^^/NDF indicates clearly that the pj exponential 
parameterization cannot properly reproduce the spectrum shape. 

For the spectra of the <p, A and A both functions give similar and acceptable ;^^/NDF. Within uncer- 
tainties, A and A have the same fit parameters. In the case of the + S spectrum, the low number 
(i.e. 3) of Pj bins cannot constrain the Levy function and therefore its ;^^/NDF in Table [5] is not de- 
fined. Nevertheless, for consistency and in order to extract particle ratios, a Levy fit is performed to 
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Table 5: Summary of the parameters extracted from the fits to the measured transverse momenta spectra using pr 
exponential ([T]l and Levy (|2]i functional forms and including point-to-point systematic uncertainties. 



Particles 


Pj exponential ([T|l 
T (MeV) ;t;2/NDF 


T (MeV) 


Levy © 
n 


;t;2/NDF 


Mesons 




325 ±4 
438±31 


117.6/14 
1.3/2 


168 ±5 
164±91 


6.6±0.3 
4.2 ±2.5 


10.8/13 
0.6/1 


Baryons 


A 
A 

+r 


392 ±6 
385 ±6 
421 ±42 


10.2/7 
5.1/7 
2.0/1 


229 ±15 
210±15 
175 ±50 


10.8 ±2.0 
9.2 ±1.4 
5.2 ±2.3 


9.6/6 
3.7/6 
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Fig. 14: Particle spectra (corrected yields) as a function of pj for Kg (circles), A (squares), A (triangles), (j) 
(stars) and + S (diamonds). The data points are scaled for visibility and plotted at the centre of the bins. 
Uncertainties corresponding to both statistics (i.e. the number of reconstructed particles) and systematics from the 
signal extraction are shown as vertical error bars. Statistical uncertainties and systematics (summarized in Table|4| 
added in quadrature are shown as brackets. The fits (dotted curves) using Levy functional form [see Eq. ©I are 
superimposed. 



obtain the integrated yields and particle ratios for all species. It must be noted that the rapidity range is 
slightly different for each species (cf. Table O. However, the rapidity dependence of particle production 
at mid-rapidity is weak enough to allow direct comparisons of the spectra [22}. 



4 Results and discussion 

The pj spectra for Kg, A, A and are shown in Fig. [14] along with the Levy fits. When comparing the 
different spectra, it is found that the inverse slope parameter T increases with the mass of the particle. 
For example, it changes from 168 ± 5 MeV for K^ to 229 ± 15 MeV for A when the Levy fit is used. 
The ± apparently do not follow this trend. However, this is most likely because the very limited 
statistics do not allow for a well-constrained fit. The shapes of the pj spectra are also compared to PHO- 
JET and PYTHIA models. For PYTHIA, several tunes (109 |201, 306 fM] and 320 r27l) are presented. 
For all species, the pj spectra are found to be slightly harder (i.e. they have a slower decrease with pj) 
than the models as presented in Figs. [151 [161 [17] and [18] For transverse momenta larger than ~ 1 GeV/c, 
the strange particle spectra are strongly underestimated by all models, by a factor of ~ 2 for Kg and even 
~ 3 for hyperons. The discrepancy is smaller in the case of the 0. 

The integrated yields {dN / dy) are obtained using the spectra in the measured range and integrating the 
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Fig. 15: Comparison of the transverse momentum ditferential yield for the Kg particles for INEL pp collisions 
with PHOJET and PYTHIA tunes 109, 306 and 320. 
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Fig. 16: Comparison of the transverse momentum differential yield for the A particles for INEL pp collisions with 
PHOJET and PYTHIA tunes 109, 306 and 320. 




Fig. 17: Comparison of the transverse momentum differential yield for the particle for INEL pp collisions with 
PHOJET and PYTHIA tunes 109, 306 and 320. 
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Fig. 18: Comparison of the transverse momentum differential yield for the S + E particle for INEL pp collisions 
with PHOJET and PYTHIA tunes 109, 306 and 320. 

Table 6: Rapidity and pj ranges, (px), corrected yields and extrapolated fraction at low pj using the Levy func- 
tion 



Particles 


l.vl 


Pj range 


(GeV/c) 


(Pt> (GeV/c) 


dN/dy 


Extrapolation (%) 


Mesons 


</• 


<0.75 
<0.60 


[0.2- 
[0.7- 


3.0] 
3.0] 


0.65 ±0.01 ±0.01 
1.00± 0.14 ±0.20 


0.184±0.002±0.006 
0.021 ±0.004 ±0.003 


12±0.4±0.5 
48± 18±7 


Baryons 


A 
A 


<0.75 
<0.75 
<0.8 


[0.6- 
[0.6- 
[0.6- 


3.5] 
3.5] 
3.0] 


0.86± 0.01 ±0.01 
0.84± 0.02 ±0.02 
0.95±0.14±0.03 


0.048 ±0.001 ±0.004 
0.047 ±0.002 ±0.005 
0.0101 ±0.0020 ±0.0009 


36±2±4 
39±3±4 
35±8±4 



Levy function for the extrapolated regions at low and high pj. The uncertainties for the dA^ /dy and (pj) 
values are computed from the errors on the fit parameters, where both the point-to-point statistical and 
systematic uncertainties of the pj spectra are taken into account. Due to the rapid decrease of the spectra, 
most of the extrapolation is done in the low pj region and amounts to 12 % for and 48 % for the (p 
(smallest and highest values respectively). Therefore, an additional uncertainty is added for the dN /dy 
to account for the uncertainty in the shape of the spectra outside the measured range: it corresponds to 
25% of the extrapolated particle yields at low pj. The measured pj ranges are specified in Table [6]for 
each particle species. 

Using the particle integrated yields presented in this paper along with the yields of charged n, K, p and 
p (16] and the measured p/p ratio [22], a comparison with STAR feed-down corrected particle ratios at 
^/s = 0.2 TeV |[25l is shown in Fig. [19] With the centre of mass energy increasing from = 0.2 TeV 
to 0.9 TeV the measured ratios are similar except the p/n^ ratio which decreases slightly from 0.068 ± 
0.011 to 0.051 ±0.005. The strange to non-strange particle ratios seem to increase but stay compatible 
within uncertainties: the K^/tt" from 0.101 ±0.012 to 0.121 ±0.013 and the A/n+ from 0.027±0.004 
to 0.032 ±0.003. 

The yields and {pj) obtained with the ALICE experiment are compared for each particle with existing 
data at the same energy and also with results at lower and higher energies. The various experiments 
differ in acceptance and event selection (i.e. NSD or INEL) but the dependence of {pj) with respect to 
these variables is found to be negligible. Consequently the {pj) values are directly comparable, whereas 
the comparison of the yields can require further scaling because of different (pseudo)rapidity coverages. 
Figure [20] reports ALICE {pj) measurements along with those of the STAR experiment |[25ll28l . It is 
remarkable that the {pj) remains close to the ISR parameterization |[29l although the collision energy 
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Fig. 19: Ratios of integrated yields including n'^^\ ¥S^\ p and p performed with the ALICE experiment llT6ll22ll 
and compared with STAR values for pp collisions at ^/i = 0.2 TeV ll25l . All ratios are feed-down corrected. For 
the ratio E^/A of ALICE, the d^V /dy|,.=o for + E is divided by 2. Statistical and systematic uncertainties are 
added in quadrature. 
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Fig. 20: {pi) vs. particle mass for the measurements performed with the ALICE experiment and compared with 
STAR values for pp collisions at ^ = 0.2 TeV Il25ll28l and the ISR parameterization ||29l . Both statistical (vertical 
error bars) and systematic (brackets) uncertainties are shown for ALICE data. 
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Table 7: The Kg mean transverse momentum and yields in INEL events from UA5, CDF, and ALICE and in 
NSD events in STAR for various ,/s. STAR results are taken from |[25 l . CDF ones and yield values with "*" are 
from ||3T1 . Other UA5 values concerning (pj) are from ll30l . 



Experiment 


V^(GeV) 


acceptance 


{pt) (GeV/c) 


dN/dy\y=o 


STAR 


200 


bl<0.5 


0.61 ±0.02 


0.134 ±0.011 


UA5 


200 


bl<2.5 


0.53 ±0.07 


0.14 ±0.02* 


UA5 


546 


bl<2.5 


0.57 ±0.03 


0.15 ±0.02* 


CDF 


630 


\y\ < 1.0 


0.5±0.1 


0.2±0.1* 


UA5 


900 


bl<2.5 


0.62 ±0.08 


0.18 ±0.02* 


ALICE 


900 


\y\ < 0.75 


0.65 ±0.01 ±0.01 


0.1 84 ± 0.002 ±0.006 


CDF 


1800 


bl < 1.0 


0.60 ±0.03 


0.26 ±0.03* 



Table 8: The (A + A) mean transverse momentum and yields for NSD events and different y/s. STAR results 
are from f25\ and UA5 results are from f32^, '331. ALICE and STAR results are feed-down corrected. The yield 
measured by ALICE has been scaled to match UA5 acceptance (\y\ < 2.0) using the method explained in section|4l 



Experiment 


V^(GeV) 


acceptance 


(Pt) (GeV/c) 


tW/dy|,v=o 


^ V+A^ per event 
measured scaled to UA5 |y| 


STAR 


200 


|y!<0.5 


0.77 ±0.04 


0.074 ±0.005 




0.24 ±0.02 


UA5 


200 


|y|<2.0 


0.80^°:?° 




0.27 ±0.07 




UA5 


546 


|y|<2.0 


0.62 ±0.08 




0.25 ±0.05 




UA5 


900 


|y|<2.0 


0.74 ±0.09 




0.38 ±0.08 




ALICE 


900 


|y| < 0.75 


0.85 ±0.01 ±0.01 


0.095 ±0.002 ±0.003 




0.46 ±0.01 ±0.02 



increased by a factor 36. Table|7]summarizes the Kg measurements performed by the UA5 Il30ll . CDF lISTTl 
and ALICE Collaborations for INEL events, and by the STAR [25 1 Collaboration for NSD events. The 
ALICE Kg yield at central rapidity, as well as the (pj), are in good agreement with UA5 results at 
900 GeV albeit with improved precision. The comparison of (A ± A) measurements are presented in 
Table [8] for NSD events. ALICE yields, measured in \y\ < 0.75 for INEL events, are scaled to the 
UA5 1321131] acceptance {\y\ < 2.0) using PYTHIA simulations. The A±A yield in NSD events is 
estimated by scaling the measured yield in inelastic events with the known ratio R of charged particle 
multiplicities in NSD and INEL events: 

/? = j^^^^^%^= 0.830 ±0.024 

This scaling factor is also used for the ALICE (E^ ± E^) yield presented in Table |9l The ALICE yields 
and (pt) for both (A ± A) and (E^ ± E^) are in good agreement with the UA5 measurements [33 1. Ta- 
ble[TO]shows the evolution of dN/dy and (pj) with the collision energy for the <p particle in NSD events. 
It includes the ALICE measurements, which are the first (p measurements at 900 GeV, and compares 
them to the results from the STAR experiment ll34l [35l at 200 GeV and the E735 experiment (3^ at 
1800 GeV. 

The baryon to meson ratio as a function of pj obtained with the (A + A) and K^ spectra measured by 
ALICE is presented in Fig. |2T] It includes the (A ± A) /2Kg ratio in pp collisions at 200 GeV measured 
by STAR 1251, and the ratios in pp collisions at 630 GeV and 1800 GeV computed with the (A + A) 
and Kg spectra pubUshed by CDF |37| and UAl l38l . UAl and CDF Collaborations provide inclusive 
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Table 9: The (E^ + S ) mean transverse momentum and yields for NSD events and different ^/s. STAR results 
are from ll25l and UA5 results are from ll33l . UA5 measures (E^ + E^) for pj > 1 GeV/c. The ALICE yield has 
been scaled to match the UA5 acceptance {\y\ < 3.0) using the method explained in section|4l 



Experiment 


(GeV) 


acceptance 


(pj) (GeV/c) 


cW/dy|,=o 


measured 


— + ) per event 
scaled to UA5 b 


STAR 

UA5 

UA5 

UA5 

ALICE 


200 
200 
546 

900 
900 


bl < 0.5 
bl < 3.0 
bl < 3.0 
bl < 3.0 
bl < 0.8 


0.90±0.01 
0.80^°i° 
L 10 ±0.02 

"■'-0.1 

0.95 ±0.14 ±0.03 


0.006 ±0.001 
0.0101 ±0.0020 ±0.0009 


0.03^°:°^ 
0.08^°:°^ 

0.05+°;°! 


0.022 ±0.006 
0.078 ±0.015 ±0.007 



Table 10: The (p mean transverse momentum and yields for NSD events and different ^/s. STAR results are 
from li34ii35J and E735 results are from li36J . The E735 Collaboration provided two values of (pj) depending on 
the functional form used to fit the data points and the uncertainties associated with each value are only statistical. 
ALICE yields measured for INEL events have been scaled to NSD as explained in section]?] 



Experiment 


^fs (GeV) 


acceptance 


{pj) (GeV/c) 


<iN/dy\y=o 


STAR 


200 


M<0.5 


0.82 ±0.03 ±0.04 


0.018 ±0.001 ±0.003 


ALICE 


900 


bl<0.6 


1.00 ±0.14 ±0.20 


0.021 ±0.004 ±0.003 


E735 


1800 


-0.4 < 3; < 1.0 


1.06±0.18 
0.94 ±0.26 


0.0186 ± 0.0041 
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CDF pp 1 .8 TeV (hl<1 .0) 
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STAR pp0.2TeV(|y|<0.5) 




UA1 pp 0.63 TeV (A+A: |ti|<2.0, K°; |r||<2.5) 



(GeV/c) 



Fig. 21: (A + A)/2K^ as a function of pj for different collision energies in pp and pp minimum bias events. The 
STAR ratio is taken from ||251 whereas the CDF and UAl ratios are computed with the (A + A) and Kg spectra 
published in ll37l and ll38l respectively. The ALICE and STAR ratios are feed-down corrected. Because the Kg 
and (A + A) spectra from UAl have incompatible binning, the Kg differential yield has been calculated for each 
(A + A) Pj data point using the fit function published by UAl. Such a choice is motivated by the fact that the X'^ 
value for the Kg spectrum fit is better than that for the (A + A) spectrum. 

spectra. The associated ratios are therefore not feed-down corrected, unlike the ALICE and STAR ones. 
The acceptance windows of these experiments differ significantly: ALICE measures A, A and in 
\y\ < 0.75, STAR in \y\ < 0.5, CDF in |t]| < 1.0, whereas UAl reconstructs (A + A) in |t]| < 2.0 and 
Kg in I T] I < 2.5. The ALICE ratio agrees very well with the STAR results in the measured pj range, 
which would suggest little or no energy dependence of (A + A) /2Kg. A similar conclusion can be drawn 
when comparing only the ratios measured by CDF at 630 GeV and 1800 GeV, although the ratio found 
by CDF for pj > 1.5 GeV/c is higher than the one observed with ALICE and STAR. The ratio computed 
from UAl spectra however shows a clear disagreement with the other measurements in an intermediate 
Pt range between pj ^ 1.5 GeV/c and pj ^ 3.0 GeV/c. PYTHIA simulations show that this discrepancy 
can not be attributed to the differences in the acceptance or in the colliding system (i.e. pp instead of pp). 



5 Conclusions 

Measurements of mesons containing strange quarks (Kg and 0) and hyperons (A, A and + S^) have 
been performed for inelastic pp collisions at ^/s = 0.9 TeV with the ALICE experiment at the LHC. The 
Levy function gives a good description of the transverse momentum spectra which have been compared 
with pQCD-based models. The Kg transverse momentum spectrum is overestimated by PYTHIA tune 
ATLAS-CSC and PHOIET below 0.75 GeV/c but is higher by a factor of ~ 2 in the pj range [1 - 
3] GeV/c. Within uncertainties, the (p meson spectrum is reasonably described by these models and the 
best agreement is obtained by PYTHIA tune D6T. We find that strange baryons are significantly under- 
predicted in both PYTHIA and PHOJET by a factor of ~ 3. The feed-down corrected ratio of baryon 
to meson as a function of pj, illustrated by the A/Kg, is consistent with the STAR measurements at 
= 0.2 TeV but lower than UAl and CDF results at = 0.63 TeV and = 1.8 TeV. The integrated 
yields and average transverse momenta have been compared with earlier data collected in pp and pp 
interactions at various energies. These results provide a useful baseline for comparisons with recent 
tunes of the PYTHIA model and a reference for future measurements in heavy-ion collisions at the LHC. 
These studies demonstrate the precision with which ALICE can measure resonances and topologically 
reconstructed weakly decaying particles. Measurements of these particles will be a substantial part of 
the ALICE programme in both pp and Pb-Pb collisions. The measurement of the (/> resonance provides 
an unprecedented reference at this energy. 
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